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We consider a possibility of identification of sources of cosmic rays (CR) of the energy above 
1 TeV via observation of degree-scale extended 7-ray emission which traces the locations of recent 
sources in the Galaxy. Such emission in the energy band above 100 GeV is produced by CR nuclei 
and electrons released by the sources and spreading into the interstellar medium. We use the data 
from the Fermi 7-ray telescope to locate the degree-scale 100 GeV 7-ray sources. We find that the 
number of such sources and their overall power match to those expected when CRs injection events 
happen every ~ 100 yr in portions of ~ 10 50 erg. We find that most of the sources are associated 
to pulsars with spin down age less than ~ 30 kyr and hence to the recent supernova explosions. 
This supports the hypothesis of supernova origin of Galactic CRs. We notice that the degree-scale 
extended emission does not surround shell-like supernova remnants without pulsars. Based on this 
observation, we argue that the presence of the pulsar is essential for the CR acceleration process. 
We expect that a significant fraction of the degree-scale sources should be detectable as extended 
sources with km 3 -scale neutrino detectors. 



I. INTRODUCTION 

The bulk of the Cosmic Ray (CR) flux hitting the 
Earth's atmosphere is composed of protons and heavy 
nuclei with energies in the 1-10 GeV range [JJ |2J. They 
are injected into the interstellar medium (ISM) by un- 
known sources and spread across the ISM via diffusion 
in the turbulent interstellar magnetic field. The diffusive 
propagation of CRs provides a major obstacle for iden- 
tification of the CR sources: angular distribution and 
spectral properties of CRs captured by particle detectors 
mostly provide information on the details of the process 
of diffusion of CRs through the Galactic magnetic field 
rather than on the location of the CR sources and mech- 
anisms of particle acceleration in the sources. 

Information on the details of diffusion of the CRs 
through the ISM, obtained from the particle detector 
measurements, provides a constraint on the energy out- 
put of the Galactic CR sources. The CR data suggest 
that typical CR particles spend only ten < 10 7 yr in the 
Galactic Disk before escaping to the larger Halo and/or 
to the intergalactic medium [HJJ. Constraint on the 
lifetime of CRs in the Galaxy imposes a constraint on 
the cumulative time-average power of the CR sources, 
Lcr ~ 10 41 erg/s [HI], needed to support the observed 
CR energy density ~ 1 eV/cm 3 [S]- 

This power could be supplied e.g. by supernovae if 
some ~ 0.1% of the energy output of the gravitational 
collapse of the star is transferred to the CRs [TJ [S] . Sev- 
eral types of the supernova- related phenomena are known 
to be associated to particle acceleration: 7-ray bursts 
(GRB), pulsars and the associated pulsar wind nebulae 
(PWN) and supernova remnant (SNR) shells. Particle 
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acceleration to the energies > 1 — 10 GeV is now firmly 
established from the 7-ray band observations of all the 
three types of objects [7J[S]. 

Alternatively, the required power could be provided by 
acceleration process which operates at the distance and 
time scales scales larger than individual supernova rem- 
nants, e.g. in superbubbles [5H11| and/or even contin- 
uously throughout the Galaxy (Ref. [12] and references 
therein) . Evidence in favor of such scenario is found from 
the study of chemical composition of the low-energy cos- 
mic ray flux [13H15] and from the direct association of 
GeV band extended emission with particle acceleration 
in Cygnus superbubble |llj . 

Clarification of the physical mechanism responsible for 
the CR production requires localization of the CR sources 
in space. This could, in principle, be done using 7-ray 
observations. This possibility arises because CRs con- 
fined in the source could interact with low-energy par- 
ticles producing 7-rays and neutrinos. However, iden- 
tification of the CR sources among the Galactic 7-ray 
sources suffers from large uncertainties. In most of the 
cases it is not clear if the amount of target material in 
the source is sufficiently large for the production of de- 
tectable CR interaction related 7-ray emission. More- 
over, in most of the sources protons and nuclei are ac- 
celerated together with electrons. The accelerated elec- 
trons produce 7-ray emission via inverse Compton and 
Bremsstrahlung mechanisms, with intensity larger than 
that of the 7-ray emission from CR interactions. As a 
result, the "hadronic" component of the 7-ray flux is hid- 
den behind much stronger "leptonic" 7-ray emission [16) . 
If CRs are injected by a process operating on large dis- 
tance scales in the Galaxy, 7-ray observations might not 
reveal individual isolated CR sources at all [12"] . 

Taking into account these difficulties, we consider in 
this paper an alternative method for the identification of 
the Galactic CR sources based on the 7-ray data. Instead 
of searching for the 7-ray emission from CR interactions 
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in the sources, we consider the 7-ray emission produced 
by CRs which escape from the source and start to spread 
into the ISM. CRs penetrating into the ISM interact with 
the ISM particles and produce 7-rays via the same mech- 
anism as inside the sources, but on much larger distance 
and time scales. The overabundance of CRs around the 
individual sources leads to the increased CR - ISM inter- 
action rate from an extended region around the source 
[TB] , This results in the appearance of extended 7-ray 
emission. We show that such extended emission around 
CR sources becomes detectable in the very-high-energy 
(VHE) band above ~ 100 GeV. We identify the extended 
Galactic sources of > 100 GeV 7-rays in the data of the 
Large Area Telescope (LAT) on board of Fermi satellite 
|17j . We study the overall distribution and energetics of 
the degree-scale sources in the Galaxy. The source statis- 
tics and morphology are consistent with the assumption 
that the extended 7-ray emission is produced by the CR 
interactions. We show that extended VHE 7-ray emis- 
sion is found around all known nearby pulsar-producing 
supernova events which occurred in the last 30 kyr, as ex- 
pected if the Galactic CRs are injected by the supernovae. 
We further argue that the presence of the young pulsars 
inside most of the degree-scale extended sources points to 
the possibility that CRs with energies above 1 TeV are 
injected by the PWN or by the composite SNR (shell- 
like SNR confining a PWN) systems, while distributed 
acceleration of cosmic rays in shocks produced by multi- 
ple supernovae and stellar wind bubbles might be most 
important at lower energies. 



where we adopt an estimate a pp ~ 4 x 10 26 cm 2 for the 
cross-section of inelastic proton-proton collisions for the 
TeV CRs PI]. 

CR interactions with ISM inside a region of the size 
r s lead to extended 7-ray (and neutrino) emission from 
the entire CR filled volume. The 7-ray luminosity can be 
readily estimated taking into account that each source 
must release some E s ~ 3 x 10 50 [IZsn /10~ 2 yr] erg 
in the form of CRs, where IZsn is the rate of CR injection 
events in the Galaxy. This leads to the estimate 
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where k ~ 0.2 is a typical fraction of the CR energy de- 
posited in the 7-rays. Note that this luminosity depends 
only on the average density of the ISM. It is largely in- 
dependent of the physical parameters of the source itself. 
The only source parameter which is fixed by the general 
energy requirement on the Galactic CR sources is the re- 
leased energy in the form of CRs, E s . The estimate of 
L 1 in Eq. Q is uncertain by at least a factor of ~ 3 
[TJ [5] . This includes the uncertainty of the estimate of 
the typical energy output per source, E s , which depends 
on the lifetime of CRs in the Galactic disk, the thick- 
ness of the disk and the rate of the injection events (with 
large uncertainties in all the three factors). Besides, the 
estimate of L 1 is energy dependent, due to the energy 
dependence of the source power and of the CR lifetime 
in the Galactic Disk. 



II. LOCATING THE CR SOURCES 

CRs which escape from sources can be traced as they 
spread into the ISM. They could be detected via the 7- 
ray emission produced in CR interactions with the ISM 
particles. CRs of the energy Ecr diffuse through the 
Galactic magnetic field. Measurements of the energy- 
dependent diffusion coefficient based on the local mea- 
surements of the primary and secondary CR nuclei and 
on the modeling of the diffuse 7-ray emission from the 
Galaxy suggest the value 
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D = D 2 s x 10 28 [E CR /4 GeV]"° cm 2 /s, 



(1) 



with an uncertainty of the prefactor D28 by up to a factor 
of ~ 3 and with the slope 8 — 0.4 ± 0.1 [3]- CRs released 
from a source some T s years ago are now filling a volume 
of size [1] 
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The average density of ISM on this distance scale is 
n iSM ~ 1 cm 3 [TB]. The typical interaction time of CRs 
propagating through the ISM is 
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III. EXTENDED VHE 7-RAY EMISSION 
AROUND CR SOURCES 

The luminosity level given in Eq. Q is high enough for 
the 7-ray emission to be detectable with existing space- 
and ground-based 7-ray telescopes. The main difficulty 
for the detection of extended emission around the indi- 
vidual CR sources is that it has to be identified on top 
of the diffuse 7-ray emission from the Galaxy [20] , 

In the energy band 0.1-1 GeV this diffuse emission ap- 
pears as a collective emission from a very large number 
of extended sources. Possible emission around individual 
sources can not be identified. The difficulty for the iden- 
tification of individual sources is, in fact, the problem 
of the source confusion. The 0.1-1 GeV 7-rays are pro- 
duced by the CRs with energies in the Ecr ~ 1 — 10 GeV 
range. In this energy range the size of the extended emis- 
sion region around the source reaches the > 100 pc scale 
in ~ 10 5 yr (see Eq. Q). Assuming the supernova rate 
of TZsn ~ 10~ 2 yr" 1 , one finds that the 100 pc scale 
regions around some N s ~ 10 3 /L> 2 8 (T^-sn/^" 2 yr -1 ) 
sources are currently filling the Galactic Disk. The den- 
sity of the sources is p s ~ N s /Vdisk ~ 3 kpc -3 , where 
V ~ 2tt(10 kpc) 2 H dlsk ~ 300 kpc 3 is the volume of the 
Galactic Disk. This implies typical distance between the 
centers of the extended emission, ~ 700 pc, so that in 
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any direction in the inner Galaxy one would find on av- 
erage several extended sources within several kiloparsec 
distance are superimposed on each other. An additional 
difficulty is the large size of the point-spread function 
(PSF) of the LAT in this energy band [T7] . 

The number of sources surrounded by the extended 
emission decreases with the increasing energy. At ener- 
gies E 1 > 100 GeV (the CR energies in the TeV range), 
the extended emission region size reaches the 100 pc in 
~ 10 kyr (see Eq. |2])). Some < 100 sources are cur- 
rently surrounded by the 100 pc scale extended emis- 
sion, so that the density of the sources decreases down 
to ~ 0.3 kpc -3 . The typical source-to-source distance is 
~ 1.5 kpc so that there is no overlap. Relatively compact 
and bright extended 7-ray emission around nearby and 
recent CR sources should start to be identifiable in the 
data. 

A source which released CRs some T s yr ago is sur- 
rounded by extended emission of the angular size 
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where R s is the distance to the source. The expected 
7-ray flux is 
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A reasonable task is to search for the relatively compact 
7-ray emission around the TeV CR sources of the age 
T s ~ 10 4 yr on top of the diffuse background. Up to 



N, 



p s TrR 2 s /2 ~ 12L> 2 



Rs 



5 kpc 



2 r 



T, 



iVa r 



10 kyr 



n 



SN 



1/100 yr 
(7) 

degree-scale extended sources could be detectable within 
the ~ 5 kpc distance in the inner Galaxy where the den- 
sity of the ISM is high enough to provide the required 
luminosity level. 

Detection / non-detection of the degree-scale sources at 
the flux level given by Eq. [6]provides a test for the super- 
nova scenario of CR production. Indeed, non-detection 
would rule out the possibility that the CRs are injected 
in portions of ~ 10 50 erg once every ~ 100 yr and would 
instead point to a scenario where the CRs are injected 
continuously throughout the Galaxy, as it is expected in 
the scenario of cosmic ray production via acceleration 
at multiple shocks formed by the stellar wind bubbles 
and supernova explosions in the star forming regions. It 
is, however, worth mentioning that in this scenario the 
degree-scale extended emission might still originate from 
the individual most active superbubbles which are cur- 
rently accelerating cosmic rays. 



IV. SEARCH FOR THE DEGREE-SCALE 
SOURCES IN THE LAT DATA 

The results obtained in the previous section suggest 
that the nearby and relatively recent sources of TeV 
CRs should reveal themselves though the degree-scale ex- 
tended 7-ray emission. The expected flux level Q is high 
enough so that the LAT sensitivity at 100 GeV is suffi- 
cient for the source detection. 

We have analyzed the data collected in the period be- 
tween August 2008 and October 2011 using the Fermi 
Science Tools 1 . The LAT event lists were filtered us- 
ing gtselect tool and only "superclean" 7-ray events 
(evtclass=4) with energies above 100 GeV arriving at 
zenith angles 9 Z < 100° were retained for the analysis. 
To identify the extended sources of the size 9 S < 1° in 
the LAT data above 100 GeV we have used the Minimal 
Spanning Tree (MST) method [2IH25] . In this method 
photons separated by angular distances smaller then a 
pre-defined value 9 max are considered as belonging to a 
" cluster" or " tree" . Point-like or extended sources could 
be found on top of the diffuse background as trees with 
large enough number of photons. Significance of the 
source detection can be estimated via comparison of the 
real data with Monte-Carlo (MC) simulations. 

To increase the sensitivity of the MST method for the 
search of extended 7-ray sources we took 9 max in the 
range 0.3° < 9 max < 0.5°, which is larger than the size 
of the PSF of the LAT telescope at the energies above 
100 GeV. In this respect the MST procedure is different 
from that used in the Ref. [23 . Close to the Galactic 
Plane, the level of diffuse sky background depends on 
the source location. The increase of 9 max leads to the 
increase of the chance coincidence probability of associa- 
tion of diffuse background photons with a tree. To esti- 
mate the background at each source location we used MC 
simulations which generate the background maps from 
the real data by changing the sky photon coordinates 
(l,b). In these MC simulations the Galactic latitude of 
the photons was left unchanged, while the longitude was 
randomly redistributed within Al = ±10° around the 
real photon position. 

We measure the extension of the sources by compar- 
ing the photon distribution around the source position 
with the PSF. To estimate the LAT PSF at E > 100 
GeV we took 1319 sources from second Fermi catalog 
[2"4] with |6| > 10° and calculated flux around them as 
function of angle. Results of this calculation are shown 
in the Fig. [I] separately for front and back photons. The 
90% containment radius is 6*90,/ = 0.14° for front pho- 
tons and #9o,& = 0.33° for back photons. We consider 
the sources for which 90% of E > 100 GeV events are 
within #90,/, 09o,& as point sources, while other sources 
are considered as extended. We remove photons associ- 
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FIG. 1: Cumulative photon distribution as function of angle 
from the source position in the energy band above 100 GeV. 
Thin lines are for PSF around extragalactic Fermi sources. 
Thick lines are photon distribution for the Galactic sources 
from Table [TTJ Solid line is for front photons, dashed line is 
for back photons. Errorbars indicate the uncertainty due to 
background fluctuations. 



ated to the point sources in the search for event clusters 
associated to extended sources. 

To identify point sources, we search for the correla- 
tion of the photon arrival directions with the positions 
of the sources from the two-year Fermi catalogue [24] . in 
a way similar to the one used in the Ref. [25] for the 
search of the extragalactic VHE 7-ray sources. The list 
of point sources found above 3er significance threshold is 
presented in the Table[TJ In this table we give the name of 
source in 2FGL catalog, galactic coordinates I and 6, the 
number of photons within 90% C.L. angle from source, 
it's significance calculated using binomial probability, the 
type of the source and its alternative name for the iden- 
tified sources. We also looked for the point sources with 
MST method with photon separation 6 = 0.2°. This 
method does not rely on the pre-defined source catalogue. 
The search with the MST method did not result in ad- 
ditional point sources which are not associated to the 
known 2FGL source. 

Ten out of 14 point sources listed in Tables [IJ are known 
VHE 7-ray sources. Four sources are new detections, 
VCS J0110+5805, 2FGL J1030.4-6015, PSR J1124-5916 
and AT20G J160350-49. Five out of 14 sources are extra- 
galactic sources found at low Galactic latitudes. Three 
sources are of unknown type. Other sources are known 
brightest Galactic 7-ray sources. All the photons associ- 
ated to the point sources listed in Table [I] are subtracted 
from the overall photon list used in the search of extended 
sources. 

The list of extended sources found using the MST 
method is given in Table [TTJ. To reduce the amount of 
false trees formed by fluctuations of the background, we 
considered only trees with known counterparts from the 



TeVCat catalogue of VHE 7-ray sources 2 for the source 
in the inner Galaxy (—60° < I < 60°). After imposing 
such a requirement, two out of 14 extended sources in 
the inner Galaxy are expected to be due to background 
fluctuations. 

Fig. [2] shows the LAT countmap in the energy range 
0.1-0.4 TeV smoothed with 0.5 degree Gaussian. The 
positions of the sources listed in Tables |T| and |TT] are shown 
on the map. 

The average radial brightness profile of the extended 
sources from Table [TTJ is shown in Fig. [T] One could see 
that the profile is clearly different from the LAT PSF. 
The cumulative photon distribution traces source pho- 
tons up to 6 ~ 0.6°. The angular sizes of the sources 
reported in the Table [TTJ are generally larger than the 
sizes reported in the HESS survey This difference 

might be attributed to the difference of the observation 
techniques used by the LAT and HESS telescopes. 



V. DETAILS FOR INDIVIDUAL EXTENDED 
GALACTIC VHE 7-RAY SOURCES 

The estimates of the section Mil show that recent 
sources of TeV CRs should be surrounded by degree-scale 
extended VHE 7-ray emission with the flux in excess of 
~ 10 -11 erg/(cm 2 s). The sources listed in Table [h] sat- 
isfy this requirement. This means that they could trace 
recent injections of CRs in the Galaxy. In this section we 
investigate the nature of CR accelerators which might 
power the observed extended emission in each individual 
source. 



A. (I, b) = (8.15, -0.14) 

This excess is associated to HESS source HESS J1804- 
216, see Fig. [3] The size of the LAT event cluster asso- 
ciated to the source is clearly larger than the extent of 
the source measured by HESS. The original discovery pa- 
per by HESS [26] suggests a PWN interpretation of the 
source associating it to a Vela-like pulsar PSR J1800-21 
[37] [55] with spin-down age T s ~ 1.7 x 10 4 yr situated at 
the distance D ~ 3 — 4 kpc [29] . 

Alternatively, an association with a shell-type SNR 
could be considered. The SNR in question is G8. 7-0.1, 
which is a relatively young (10 yr) and nearby SNR vis- 
ible in the 0.1-100 GeV energy band [30 . Possible asso- 
ciation of G8.7-0.1 with PSR J1800-21 is considered (see 
e.g. [27]) based on the proximity of the two objects in 
space and similar age. 
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FIG. 2: LAT countmap of the inner Galaxy in the 100-400 GeV band smoothed with a 0.5 degree Gaussian. VHE 7-ray sources 
from Table III] are shown by green circles. 



B. 



(l,b) = (17.58,-0.14) and (1,6) = (16.74,-1.79) 



This extended emission region includes the HESS 
source HESS J1825-137, which is identified with the pul- 
sar PSR B1823-13 [311 132], see Fig. E Similarly to PSR 
J1800-21, this is a relatively nearby (4 kpc) and young 
(20 kyr) pulsar [59] . 

The extended radio source, possibly a supernova rem- 
nant, G 17.4-0.1 [35], is found within the source exten- 
sion, so that relation of the extended emission to a shell- 
like SNR could not be ruled out. Lack of basic informa- 
tion on G 17.4-0.1 (age, distance) makes the analysis of 
association of the observed extended emission to G 17.4- 
0.1 impossible. 

Immediately adjacent to the extended source around 
HESS J1825-137, there is an extended emission which 
includes a 7-ray loud binary LS 5039 (Fig. [4|. It is 
not clear if this source is just a further extension of a 



brighter HESS J1825-137, or it is an independent source. 
The relation of this source to LS 5039 is also not evident. 



Energy-dependent morphology of HESS J1825-137 was 
studied in detail in the Ref. [35]. Detailed observations 
by HESS have revealed a more compact emission region 
of the size ~ 0.3° at the highest energies on top of a 
more extended degree-scale emission at the energies of 
several hundred GeV. Suzaku observaitons of the HESS 
J1825-137 region reveal an extended X-ray emission asso- 
ciated to the central more compact 0.3° source, but not 
to the extended degree-scale source [31] . The X-ray emis- 
sion detected by Suzaku could be consistently interpreted 
as synchrotron emission from electrons in the PWN of 
PSR B1823-13. The absence of X-ray emission from the 
degree-scale source might indicate that the degree-scale 
7-ray emission is due to CR interactions. 
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TABLE I: List of point like VHE sources found at the positions of known Fermi sources. I, b are source coordinates, Nph is the 
number of photons within the 90% containment radius of the PSF, P is the chance coincidence probability for the event cluster 
to be a background fluctuation. Notations in the column "Type" are PWN: pulsar wind nebula; SNR+CCO: supernova remnant 
with compact central object, BLZ: blazar; AGU: AGN of unknown type; GRLB: 7-ray loud binary. Previously reported VHE 
7-ray sources are marked in bold. Other sources are new detections in the VHE band. 




FIG. 3: Sky region around extended source at (l,b) = 
8.15,-0.14). White circle marks the position of 2FGL ex- 
tended source associated to supernova remnant W30. 



FIG. 4: Sky region around extended source at (/,&) = 
17.58,-0.14). Red circles show positions and sizes of SNRs 
from the Green catalogue [42] , 



C. (l,b) = (23.32,-0.16), (25.21,-0.16) and 
(26.87,-0.12). 

The event clusters associated to these three adjacent 
excesses merge into one extended emission region along 
the Galactic Plane, in which the central source, the clus- 
ter at (l,b) — (25.21 — 0.16), dominates. This central 
source is associated to the HESS source HESS J1837-069 



identified with a PWN of the pulsar PSR J1838-0655 
[331 13"6"] , which has the spin-down age of ~ 22 kyr [2H] ■ 

The side sources are also associated to the HESS 
sources HESS J1834-087 and HESS J1841-055 (Fig. |). 
The most probable low-energy counterpart of HESS 
J1834-087 is the SNR W41, which contains a central 
PWN-like object CXOU J183434.9-084443 [37||3S]. Ref. 
[3"9"] provides an estimate of the age of this SNR, T s ~ 
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0.23 


0.37 


l.e-6 


11 


4.6 ± 1.4 


HESS J1857+026 




J1856±0245 


10.3 


2.0 


8 


78.09 


2.54 


0.33 


0.38 


l.e-5 


7 


2.3 ±0.9 


VER J2019+407 


7Cyg 


J2021±4026 




7.7 


9 


284.32 


-0.57 


0.32 


0.42 


7.e-3 


4 


1.3 ±0.7 


Westerlund 2 




J1023-5746 




0.5 


10 


287.12 


-0.80 


0.46 


0.74 


2.e-4 


9 


2.9 ± 1.0 


near Eta Car 










11 


313.56 


0.11 


0.2 


0.32 


8.e-6 


8 


2.6 ± 1.0 


Kookaburra 




J 1420-6048 


7.7 


1.3 


12 


331.66 


-0.58 


0.27 


0.64 


7.e-4 


11 


3.7 ± 1.1 


HESS 1614-518 




J1614-5144 






13 


332.57 


-0.18 


0.34 


0.63 


l.e-3 


10 


3.3 ±1.0 


HESS J1616-508 




J1617-5055 


6.5 


0.8 


14 


336.25 


0.04 


0.37 


0.59 


l.e-6 


16 


5.4 ± 1.3 


HESS J1632-478 




J1632-4757 


7.0 


24 


15 


339.56 


-0.79 


0.37 


0.72 


3.e-3 


10 


3.4 ± 1.0 


Westerlund 1 




J1648-4611 


5.7 


11 


16 


344.90 


0.23 


0.72 


1.05 


3.e-2 


8 


2.8 ± 1.1 


HESS J1702-420 




J1702-4128? 


5.2 


5.5 


17 


346.20 


-0.31 


0.37 


0.57 


l.e-2 


7 


2.7 ± 1.0 


HESS 1708-410 




J 1706-4009? 


3.8 


0.9 


18 


358.06 


-0.54 


0.57 


0.63 


l.e-4 


10 


3.7 ± 1.2 


HESS J1745-303 











TABLE II: List of extended sources found in the MST analysis. #50, #90 are the radii within which 50% and 90% of the 
tree photons are contained. N p h is the background subtracted number of photons in the tree. P90 is the chance coincidence 
probability of background fluctuation which would have N p h photons within a circle of the radius #90. F is the source flux 
in units of 10 -11 erg/cm 2 s. SNR and PSR columns give the names of supernova remnants and/or pulsars suggested to be 
associated to the source. R s and T s are the source age (in the units of 10 4 yr) and distance (in kpc). 
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FIG. 5: Sky region around extended sources at (I, b) = 
(23.32,-0.16), (25.21,-0.16) and (26.87,-0.12). Green 
cross inside HESS J1834-087 marks the position of CXOU 
J183434. 9-084443. Red circles show positions and sizes of 
SNRs from the Green catalogue |42) . 



10 kyr. The extension of the cluster of events associated 
to HESS J1834-087 is, however, larger than the size of 
the HESS source. Within this extension another rela- 
tively young pulsar is present, PSR B1830-08 (see Fig|5]), 
which somewhat younger, T s ~ 30 kyr. 

The source associated to the HESS J1841-055 does not 
have a shell-like SNR candidates within its extension. At 
the same time, a young nearby pulsar PSR J1841-0524 of 
spin-down age T s ~ 30 kyr |29j is situated in the center of 
the source. This pulsar could be considered as a possible 
low energy counterpart of the extended 7-ray source. 



D. (I, b) = (36.20,0.02) 

This extended source is spatially coincident with HESS 
source HESS J1857±026 and is identified with a PWN of 
PSR J1856+0245 [Ml IS] of spin-down age T s ~ 20 kyr. 
No shell-like supernova remnant from the Green cata- 
logue [12] is found within the source extent (see Fig. |6j). 
The size of event cluster associated with the source is 
much larger than the reported size of the HESS source 
(see Fig. [61). The source is elongated in the direction of 
another HESS source, HESS J1858±020. 




FIG. 6: Sky region around extended source at (I, b) = 
(36.20,0.02). Red circle shows position and size of a SNR 
from the Green catalogue |42| . 




FIG. 7: Sky region around extended source at (l,b) 
(78.09,2.54). Red circle shows the extent of iCygni SNR. 



E. (I, b) = (78.09, 2.54) 

This source is a part of extended emission from the 
PSR J2021+4026 - 7Cygni PWN+SNR system. The 
centroid of extended emission is displaced from the pul- 
sar position, but the emission still comes from within the 
SNR boundary (Fig. [7|), so that it is not the extended 
emission from the CRs escaping into the ISM. Taking 
into account the age of the SNR, T s ~ 80 kyr, and its 



proximity R s ~ 1.5 kpc [43 one could find that the ex- 
pected size of the extended emission from the escaping 
CRs is 6 S ~ 5° ([5]), i.e. it is expected to cover the entire 
Cygnus region. 

GeV band extended emission on the angular scale ~ 5° 
from the sky region immediately adjacent to the 7Cygni 
source was recently discussed in the Ref. [11] . This refer- 
ence attributes the extended emission to the interactions 
of CRs in the ISM. Although the possibility of relation 
of this extended emission to 7Cygni is mentioned in this 
reference, an alternative possibility for the origin of CRs 
interacting in the ISM is considered, based on the spa- 
tial displacement of the source from the 7Cygni position. 
The scenario favored in the Ref. |llj is injection of cosmic 
rays via acceleration at multiple shocks formed by stel- 
lar wind bubbles in the OB associations in the Cygnus 
region, which form the Cygnus superbubble. 

We notice that in principle the displacement of the 
extended emission from the possible primary source po- 
sition (pulsar) is expected, because of the presence of 
the ordered component of Galactic magnetic field. Or- 
dered magnetic field introduces anisotropy in the CR dif- 
fusion from the source, so that the extended emission pro- 
duced by CRs escaping from the source is not expected 
to be symmetrically distributed around the source po- 
sition. Additional asymmetry could be introduced by 
spatial variations of the ISM density. 

It is also worth noticing that, contrary to the 
low-energy (GeV) emission, the maximum of the ex- 
tended emission in the VHE band is spatially coincident 
with the 7Cygni supernova remnant associated to PSR 
J2021+4026. This suggests that the higher energy (TeV) 
and lower energy (1-10 GeV) might be produced at dif- 
ferent locations, via different acceleration mechanisms. 



F. (I, b) = (284.32, -0.57) and near Eta Car 

The excess at (I, b) = (284.32 - 0.57) is positionally 
coincident with the HESS source HESS 1023-575 identi- 
fied as a PWN of PSR J1023-5746 (previously associated 
with Westerlund 2 cluster) [33]. The pulsar is relatively 
young, T s ~ 5 kyr. Absence of the radio emission does 
not allow the measurement of the distance based on the 
dispersion measure 29J . No shell- like supernova remnant 
from the Green catalogue is visible within the source ex- 
tension. 

No obvious counterpart candidates could be found 
within the extent of the emission region adjacent to the 
Eta Car (see Fig. [8]): no known pulsars or shell-like 
SNRs. The nature of this excess requires further inves- 
tigation. The emission from this sky region comes from 
the direction tangential to the Carina-Sagittarius arm of 
the Galaxy (see Fig. 14 below). The enhancement of the 
100 GeV emission might, in principle, be an enhancement 
of Galactic diffuse emission due to the projection effect. 
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FIG. 10: Sky region around extended source around 
(l,b) = (331.66,-0.58), (l,b) = (332.57,-0.18), (l,b) = 
(334.97,-0.26) and (l,b) = (336.56,0.04). Red circles show 
positions and sizes of SNRs from the Green catalogue 



FIG. 8: Sky region around extended source at (l,b) == 
(284.32, —0.57) and near Eta Car. Red circle shows position 
and size of a SNR from the Green catalogue [42| . 



H. (I, b) = (331.66, -0.58), (I, b) = (332.57, -0.18) and 
(l,b) = (336.56,0.04) 




FIG. 9: Sky region around extended source around (l,b) 
(313.56,0.11). 



G. (I, b) = (313.56,0.11) 



The extended source at (l,b) = (313.56,0.11) is spa- 
tially coincident with the known HESS source at the 
Kookaburra PWN around PSR J1420-6048 of the age 
13 kyr situated at the distance 7.7 kpc [29] . 



The three excesses are close to each other so that the 
event clusters associated to these three sources merge 
into a single extended emission region along the Galac- 
tic Plane. The excesses could be identified with known 
HESS sources: HESS J1614-518, HESS J1616-508 and 
HESS J1632-4757 (see Fig. 
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The source at the position of HESS J1614-518 does 
not have any known pulsars or SNRs within its extent. 
The direction toward this source is tangent to the Norma 
arm of the Galaxy. Similarly to the case of the Scrutum- 
Centaurus arm, enhancement of the matter density in 
the arm and projection effects might be responsible for 
the enhanced VHE 7-ray flux from this direction. 

HESS J1616-508 is most probably associated to the 
PWN of PSR J1617-5055 g5] of the age T s ~ 8 kyr 
[29], although several shell-like SNRs are adjacent to the 
HESS source position and are within the extent of the 
LAT event cluster at (l,b) = (332.57,-0.18) (see Fig. 
101. 



The event cluster at the position of HESS J 1632-478 is 
the brightest excess in this source group. The maximum 
of excess is at the position of the pulsar PSR J1632-4757 
and Ref. [46] associates the source to the pulsar. In this 
case the age of the supernova associated to the source is 
T s ~ 24 kyr. Several other pulsars and SNRs, as well as 
another HESS source, HESS J1634-472, are found at the 
periphery of the source (see Fig. 10). It is possible that 



the overall emission from this bright excess is composed 
of contributions from several sources. 
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, HESS, J17G2-J 



RXJ1713.7-394* HESSV1T08-410 U , 



FIG. 11: Sky region around extended source at (l,b) 
(339.56,-0.79). 



FIG. 12: Sky region around extended sources at (I, b) = 
(344.90,0.23) and (l,b) = (346.20,-0.31). Red circles show 
positions and sizes of SNRs from the Green catalogue [42] . 



I. Source at (l,b) = (339.56,-0.79) 

Excess at this position coincides with the HESS source 
in the direction of Westerlund 1 stellar cluster. The LAT 
event cluster associated to the source is more extended 
than the HESS source toward the position of the LAT 
pulsar PSR J1648-4611 (Fig. pj. This is a T s ~ 100 kyr 



pulsar situated in the Norma arm of the Galaxy at the 
distance R s ~ 5 kpc. 



J. (I, b) = (344.90, 0.23) and (I, b) = (346.20, -0.31) 

The LAT event cluster in this direction extends along 
the Galactic Plane between three HESS sources, HESS 
J1702-420, HESS J1708-410 (Fig. 



12 1. Several shell-like 



SNR and PWNe are found within the source extent, so 
that no "natural" counterpart candidates could be iden- 
tified. Maximum of the excess at 100 GeV is at HESS 
J1708-410, but at slightly lower energies the maximum 
shifts toward HESS J1702-420, close to the pulsar PSR 
J1702-4128, so that the source morphology is energy- 
dependent. This pulsar has the age 55 kyr and is situated 
at the distance ~ 5 kpc [29] . Another pulsar, PSR J1706- 
4009, is situated within the extent of the event cluster at 
(l,b) = (346.20,-0.31). This pulsar is situated at the 
distance 3.8 kpc and has the spin-down age 9 kyr 29 . 




FIG. 13: Sky region around extended source at (l,b) = 
(358.06, —0.54). Red circles show positions and sizes of SNRs 
from the Green catalogue 42 . 



within the source extension so that no obvious counter- 
part candidate could be singled out. 



K. Source at (I, b) = (358.06, -0.54) 

This excess is at the position of the HESS source HESS 
J1745-303, but has a much larger size than the HESS 
source. Several pulsars and shell-type SNRs are found 



VI. THE NATURE OF EXTENDED SOURCES 



The overall number of extended 100 GeV 7-ray sources 
at the flux level in excess of 10 -11 erg/cm 2 s is consis- 
tent with expected number of the regions of degree-scale 
extended emission around recent TeV CR sources in see- 
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FIG. 14: Locations of the pulsars associated to the extended 
sources at 100 GeV on the Galaxy map. Boxes mark the 
brighter excesses identified in the Fermi data, blue circles 
correspond to pulsars within weaker VHE 7-ray sources de- 
tected by the ground-based 7-ray telescopes. Red crosses 
mark the positions of known pulsars with ages in the range 
T a < 3 x 10 4 yr. Blue circle marks the position of the Galac- 
tic Center. Blue box is the position of the Sun. Solid curves 
show the locations of the main arms of the Galaxy. Dashed 
arcs show the distances 5 kpc and 8 kpc. 



nario in which CRs are injected in portions of ~ 10 50 erg 
rather than continuously all over the Galaxy (see sec- 
tions [n III ) . The source sizes are also consistent with an 
assumption that the extended emission is produced by 
high-energy particles diffusing into the ISM. Fifteen out 
of the eighteen extended sources listed in the Table |TT] 
are possibly or clearly associated with pulsars of the age 
T s < 30 kyr. This is consistent with the hypothesis that 
TeV CRs are injected by supernova related phenomena. 
Indeed, majority of the core collapse supernova (more 
than 75%) in the Galaxy are expected to result in pro- 
duction of the neutron star [47] . Most of the collapsing 
massive stars are rotating, so that significant fraction of 
the new-born neutron stars are expected to be visible as 
pulsars [4*5] . 

An essential element of the supernova scenario for the 
origin of CRs is that a significant part of supernovae 
(rather than just a small fraction of them) should inject 
some 10 50 erg of CRs in the interstellar medium. As- 



sumption that every supernova injects ~ 10 erg in CRs 
implies that the 100 pc scale extended VHE 7-ray emis- 
sion regions should exist around nearly all the relatively 
recent (30 kyr) supernovae. Calculations of the previous 
sections suggest that extended emission regions situated 
at the distance R s < 5 kpc should be detectable at the 



flux level in excess of 10 -11 erg/cm 2 s. 

The sensitivity of LAT in the E > 100 GeV energy 
range is only marginally sufficient for detection of degree- 
size sources with fluxes ~ 10~ n erg/cm 2 s. In fact, all 
the sources listed in Table |H] have slightly higher fluxes. 
This means that, in principle, some extended emission 
regions around the nearby TeV CR sources of the age T < 
30 kyr might have escaped detection by LAT. At the same 
time, their flux level might be above the sensitivity of the 
ground-based 7-ray telescopes. To verify this conjecture, 
we compiled a list of all known pulsars with the spin- 
down age shorter than 30 kyr using the ATNF pulsar 
catalogue [2S]. As it is explained above, such pulsars 
could be used as tracers of recent supernova explosions. 
Fig. [l4] shows the locations of the pulsars in the Plane of 
the Galaxy. There are 56 such pulsars at the distances 
smaller than 20 kpc. Some 29 of them are found within 
0.5 degrees from the known VHE 7-ray sources. The 
list the pulsars found in vicinity of known VHE sources, 
which are not detected in the MST analysis of LAT data 
is given in Table |III| In fact, a large fraction of the sources 
listed in Table |III| are associated with 100 GeV 7-ray 
event clusters in LAT, but with the event statistics below 
the threshold adopted in our analysis. 

The positions of the pulsars associated with the bright- 
est extended VHE 7-ray sources listed in Table [H] are 
marked by magenta boxes in Fig. |14| One could notice 
that most of the sources are situated in the Norma arm 
of the Galaxy, which is ~ 5 kpc away from the Sun. The 
fact that the extended sources in the Norma arm appear 
brighter than the extended emission around closer CR 
sources might look surprising at the first sight. How- 
ever, the Norma arm is a special region of the Galaxy 
with the highest star formation rate [49]. High star for- 
mation rate is associated to the increased matter den- 
sity in the arm (presence of numerous molecular clouds) , 
which boosts the VHE 7-ray luminosity of the extended 
emission. Positions of the pulsar-associated extended or 
point sources which are below the LAT sensitivity above 
100 GeV, but possibly associated to the VHE 7-ray emis- 



sion are marked by the blue circles in Fig. 14 
From Fig 



14 one can see that most of the known 



nearby supernovae associated with pulsars of the age 
T s < 30 kyr are spatially coincident with excess VHE 7- 
ray emission. In the case of sources listed in Table [TTJ the 
emission is extended, with the typical source size in the 
degree range. This implies the spatial extent ~ 100 pc 
at the source distance. Table Mil includes both isolated 
(more compact) sources associated to the young pulsar 
wind nebulae (like e.g. MSH 15-52) and extended emis- 
sion regions discussed above. Further investigation of the 
nature of emission from these sources is needed. 



VII. DISCUSSION 

Existence of the 100 pc size extended VHE 7-ray emis- 
sion regions around the locations of recent injections of 
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Name 


PSR 


R s 


T 

-L S 


1 


VelaX 


B0833-45 


0.29 


1.1 


2 


G292.2-0.5 


J1119-6127 


8.40 


0.2 


3 


HESS J1303-631 


J1301-6305 


15.84 


1.1 


4 


HESS J1356-645 


J1357-6429, 


4.09 


0.7 


5 


Rabbit 


J1418-6058 




1.0 


6 


MSH 15-52 


B1509-58 


5.8 


0.2 


7 


HESS J1708-443 


B1706-44 


1.82 


1.8 


8 


HESS J1741-302 


B1737-30 


3.28, 


2.1 


9 


GO.9+0.1 


J1747-2809 


> 8 


0.5 


10 


HESS J1809-193 


J1809-1943 


3.57 


1.1 






J1811-1925 




2.3 


11 


HESS J1813-178 


J1813-1749 




0.5 


12 


HESS J1833-105 


J1833-1034 


4.30 


0.5 


13 


HESS J1846-029 


J 1846-0258 


5.10 


0.1 


14 


MGRO J1908+06 


J1907+0602 


3.01 


2.0 


15 


G54.1+0.3 


J1930+1852 


5.00 


0.3 


16 


MGRO J2019+37 


J2021+3651 


18.9 


1.7 


17 


Boomerang 


J2229+6114 


3.0 


1.1 



TABLE III: List of T s < 30 kyr pulsars positionally coincident 
with known VHE 7-ray sources 



TeV CRs in the Galaxy is an immediate testable predic- 
tion of the scenario of injection of CRs by the supernovae. 
These extended emission regions should be visible as de- 
gree scale VHE 7-ray sources with fluxes at the level of 
~ 10~ n erg/cm 2 s above 100 GeV. In this paper we have 
searched for such degree-scale sources in the LAT data. 
We have found 18 such sources. Table |H] gives the list of 
the extended sources found using the MST method. The 
overall extended source statistics, as well as the morpho- 
logical properties of sources listed in Table [TT] agree with 
the theoretical expectations. Based on this agreement, 
we suggest that the bright extended sources are associ- 
ated to the sites of recent CR injections in the Galaxy. 

The largest part (80%) of the extended VHE 7-ray 
sources in the Galactic Plane, listed in Table |IT| could be 
associated to pulsars. For eleven out of fifteen pulsar- 
related sources (i.e. 70%) no shell-like SNR structure 
has been found, which might signify that the supernovae 
did not result in production of a regular expanding shell 
(similarly e.g. to the supernova which produced the Crab 
pulsar). Thus, it is not clear if the presence of a regu- 
lar shell is an essential ingredient of the CR acceleration 
process. At the same time, all the shell-type SNRs enter- 
ing Table [TT] are composite SNRs confining a PWN. None 
of the sources are associated to shell-like SNR without 
a pulsar. This might indicate that CR acceleration phe- 
nomenon is related to the presence of a pulsar in the 
SNR. 

A significant fraction of pulsars is born in OB associ- 
ations [50l [51] . Young pulsars which have typical birth 
kick velocities < 10 3 km/s= 1 pc/kyr do not leave their 
parent association of the size ~ 10 — 100 pc until they 



reach the age 10 4 — 10 5 yr. Taking this into account, 
it is possible that locations of the young pulsar in the 
Galactic Plane just trace the OB associations where they 
were born. In this case the degree-scale extended 7-ray 
emission might, in fact, be related to the presence of an 
OB association or a superbubble, which contains multi- 
ple supernova remnants and bubbles blown by the winds 
of massive stars. The most clear example of this type 
is given by the extended source in the Cygnus region 
(source number 8 in Table]!!]) . Study of extended emission 
in the GeV band reported Ref. [TT] finds that cosmic rays 
were recently injected in this region and suggests that the 
injection mechanism is related to the multiple shocks pro- 
duced by the collective effect of winds from massive star 
filling the Cygnus region superbubble. An alternative 
possibility, also discussed in the Ref. [TT] is that the cos- 
mic rays are produced by the nearby 7Cygni supernova 
which hosts PSR J2021+4026 (see section |VE|. 

Spatial coincidence of the extended VHE 7-ray sources 
in the Galactic Plane with the pulsars introduces a po- 
tential possibility that the extended emission is produced 
by electrons escaping from the PWN around the pul- 
sar, rather than by CRs filling the 100 pc-scale bub- 
ble. This possibility requires further investigation. The 
VHE 7-ray emission in the E ~ 100 — 300 GeV range 
could be produced via inverse Compton scattering by 
electrons of comparable energies upscattering optical- 
infrared Galactic photon background (photons with en- 
ergies ~ 0.1 — 1 eV). The inverse Compton cooling time 
of electrons is ~ 2 x 10 6 yr, assuming the density of the 
optical-infrared background U ~ 1 eV/cm 3 , typical for 
the Norma Arm [S5]. This cooling time scale is order- 
of-magnitude shorter than the CR cooling time pi), cal- 
culated assuming the ISM density n ~ 1 cm" . This 
implies that electrons loose energy more efficiently than 
protons. To compare the luminosities generated by elec- 
trons to those produced by CRs, one needs to know the 
CR and electron densities in the relevant energy scales, 
which is highly uncertain. A first guess could be found 
from the ratio of electron and CR fluxes observed at the 
Earth: the electron flux at 300 GeV is - 10~ 2 of the TeV 
CR flux [53, 54 j. Assuming this ratio, one finds that the 
CR contribution to the extended VHE source flux should 
dominate over the electron contribution. It is, however, 
not clear if the locally observed electron-to-nuclei ratio 
could serve as a good estimate of the ratio close to the dis- 
tant CR sources. Electrons should still give a dominant 
contribution to emission from the central parts of the de- 
gree scale sources, at the locations of Crab-like PWNe 
around the pulsars. Two complementary measurements 
could be used to separate electron and CR contribution 
to the VHE 7-ray flux from the degree-scale sources. 

First, electron-dominated emission from the central 
PWN could be identified via radio-to-X-ray observa- 
tions which could establish the morphology of the PWN 
and measure the spectrum of the synchrotron emis- 
sion from high-energy electrons. At larger, degree 
scales, electrons with energies E e ~ 300 GeV are ex- 
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pected to produce synchrotron emission in the e s ~ 
1(T 2 [B/3 fiG] [£ e /300 GeV] 2 eV energy range, i.e. in 
the far infrared (assuming magnetic field B ~ 3 /iG, typi- 
cal for the ISM) . Detection of the infrared counterparts of 
the extended VHE 7-ray sources will constrain the elec- 
tron contribution to the degree-scale VHE source flux. 

Alternatively, proton/nuclei contribution to the ex- 
tended source flux could be measured via detection 
of neutrino signal from the sources listed in Table [TT] 
with km 3 scale neutrino detectors like IceCube [55] and 
km3net [56] . Asuming that most of the 7-ray emission 
from sources listed in Table [TT] is produced by CR pro- 
tons/nuclei, the neutrino flux is expected to be compara- 
ble to the 7-ray flux given in the column 8 of the Table. 
This would make the sources from Table [IT] the strongest 
Galactic neutrino sources. The slope of the neutrino 
spectrum could be inferred from the HESS observations 
|26j . The neutrino spectrum produced in nuclei- nuclei 
interactions is expected to follow the 7-ray spectrum, i.e. 
to have the powerlaw shape dN v /dE ~ E~ 2A in the en- 
ergy range up to at least ten(s) of TeV. This slope is 
harder than the slope of the atmospheric neutrino spec- 
trum, so that the sources should be preferably detectable 
at the highest energies. However, one should note that 
the source size is expected to increase with energy (see 
Eq. 5) making the sources extended for IceCube and 
km3net. 

IceCube is mostly sensitive to the sources in the North- 
ern hemisphere. Only two of the 18 sources in Table [IT] 
are situated in the Northern hemisphere: the source 8 in 
the Cygnus region and the source 9 associated to HESS 
J1857+026. Currently existing upper limits on the point 
source fluxes with IceCube 40-string configuration [57] at 
positions of these sources are still above the expected neu- 
trino flux level, expected to be around ~ 10~ n erg/cm 2 s 
at the energy of 0.1 TeV and still lower at the energies 
above 10 TeV at which IceCube is most sensitive. 

Possible association of the degree-scale VHE 7-ray 
sources to pulsar/PWN systems, rather than to the shell- 
like SNRs without pulsars suggests a possibility that the 
presence of pulsar plays a significant role in the CR 
production. PWNe formed by the collisions of pulsar 
winds with the SNR material and/or with the interstellar 
medium are known to be very efficient particle accelera- 
tors. In these systems most of the energy initially stored 
as the rotation energy of the neutron star is transferred to 
the high-energy particles. Most of our knowledge of par- 
ticle acceleration mechanisms in PWNe is derived from 
the observations of radio-to-7-ray emission from the ac- 
celerated electrons. In general, the standing relativistic 
shock at the interface of the pulsar wind and the super- 
nova/interstellar medium should also accelerate protons 
and heavier nuclei. The energy stored in the rotation of 
the neutron star is 
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for a typical neutron star with the moment of inertia 



/ ~ 1.5 x 10 45 g cm 2 and the angular velocity = 
2ir /Pini corresponding to the initial rotation period Pini- 
If typical initial periods of the pulsars at birth are in the 
range Pj„j < 10 ms, the rotation energy of the neutron 
star is sufficient for powering the CR production. 

Relativistic nature of the shock accelerating CRs could 
leave an imprint on the CR spectrum in the form of a 
hardening below the energy corresponding to the gamma 
factor of the shock [58 . Such a hardening might be re- 
sponsible for the low-energy break observed in the spec- 
trum of Galactic CRs [5]. 

Periods of known pulsars increase with time as a pow- 
erlaw P ~ where n is the breaking index which, 
for most of the pulsars is in the range 2 < n < 3 [59] . 
Observations of the powerlaw spin-down do not provide 
information on the initial rotation period of the neutron 
star at birth. The initial periods of rotation of neutron 
stars have to be estimated based on numerical model- 
ing of stellar evolution, supernova explosions and of the 
subsequent evolution of the neutron stars. Numerical cal- 
culations of evolution of massive stars up to the point of 
gravitational collapse only now start to be able to ac- 
count for the effects of the stellar rotation. The state- 
of-art simulations predict the initial period in the range 
of Pi n i ~ 10 ms for a typical star of 15 Solar masses, if 
a special mechanism for the removal of the angular mo- 
mentum from the stellar core by strong magnetic field 
is assumed 48J. Otherwise, the calculations suggest the 
initial periods in the range of Pini ~ 1 ms in the absence 
of strong magnetic coupling of the core to the outer layers 

This implies that the rotation energy of the neutron 
star could, in principle, provide the power sufficient for 
production of CRs. If this case the presence or absence of 
a well-defined SNR shell is not of crucial importance for 
the acceleration. If the shell is present, the high-energy 
particles accelerated in the PWN would be confined by 
the magnetic fields in the shell, where they could interact, 
like it happens in the MSH 15-52 / PSR B1509-58 sys- 
tem. Otherwise, if the regular shell is not formed (like in 
the Crab nebula), particles would escape from the PWN 
directly into the interstellar medium. Alternatively, pul- 
sars might play an " auxiliary" role in the CR acceleration 
process providing an initial injection of "pre-accelerated" 
particles in the SNR shell. Understanding of the role of 
pulsars in the CR acceleration process requires further 
investigation. 

Analysis of abundances of elements in the cosmic ray 
flux suggests that the bulk of Galactic cosmic rays with 
energies below 10 GeV is accelerated in the environments 
containing a mix of supernova ejecta and older interstel- 
lar material typically present in OB associations and su- 
perbubbles [J3]rri]. Moreover, the the time scale of accel- 
eration is most probably longer than 100 kyr, an obser- 
vation based on the observation of abundances of specific 
elements [15]. This implies that low-energy cosmic rays 
are not produced out of the single supernova ejecta soon 
after the supernova explosion, but rather originate from 
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the slow acceleration on the system of multiple shocks 
present in OB associations and superbubbles. Recent 
Fermi observations of 7-ray emission from the interac- 
tions of freshly accelerated cosmic rays in the Cygnus 
superbubblc further support this scenario. 

At the same time, recent observations of several breaks 
in the interstellar cosmic ray spectrum: softening above 
- 10 GeV [5J and hardening above 200 GeV [53] indi- 
cate the presence of multiple components in the cosmic 
ray spectrum. The component which dominates the low- 
energy cosmic ray flux and peaks in the 1-10 GeV range 
has soft spectrum above 10 GeV and becomes subdom- 
inant above 200 GeV. It is not clear if conclusion about 
possible OB associations and/or superbubble origin of 
the dominant low-energy component holds also at the 
energies above 200 GeV, or a different acceleration mech- 
anism is relevant at higher energies. As it is mentioned 
above, association of the degree-scale extended sources 
with young pulsar might reflect possible OB association 
membership to the pulsar parent stars, so that multiple 
acceleration mechanisms related to the pulsar wind neb- 
ulae, supernova remnants and to the entire star-forming 
region around the supernova/pulsar might be responsible 
for the production of cosmic rays filling the 100 pc-size 
regions visible as the desgree-scale 7-ray sources. 

VIII. CONCLUSIONS. 

We have considered a possibility of localization of the 
sources of TeV CRs via detection of degree-scale ex- 
tended VHE 7-ray emission produced by CRs escaping 
into the ISM. We have identified extended sources of 



E > 100 GeV 7-rays in the Galactic Plane using the 
data of LAT and studied their statistics and morpho- 
logical properties. The list of the degree-scale extended 
source is given in Table [TT] The average surface brightness 
profile of the extended sources is shown in Fig. [I] 

The observed extended source properties suggest their 
interpretation as extended emission produced by CRs es- 
caping from relatively young sources of TeV CRs. Most 
of the degree-scale extended sources are associated to 
nearby pulsars with the spin-down age T < 30 kyr, as 
expected in the model of supernova origin of the TeV 
CRs (if pulsars are considered as tracers of the recent su- 
pernova events) and their star-forming regions. We 
find that nearly all known young nearby pulsars are as- 
sociated to excess extended VHE 7-ray emission. At the 
same time, none of the detected degree-scale VHE sources 
is associated to a SNR without a pulsar. 

The observed degree-scale extended emission could be 
produced by both the proton/nuclear CRs and elec- 
trons/positrons escaping from the sources. It is difficult 
to distinguish between the two contributions based on the 
7-ray data only. We suggest that the two contributions 
could be distinguished based on the complementary ob- 
servations in the infrared band and in the VHE neutrino 
channel. 
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